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SUMMARY 


A small— deflection theory is developed for the elastic "behavior of 
orthotropic flat plates in which deflections due to shear are taken into 
account. In this theory, which covers all types of flat sandwich 
construction, a plate is characterized by seven physical constants (five 
stiffnesses and two Poisson ratios) of which six are independent. Both 
the energy expression and the differential equations are developed. 
Boundary conditions corresponding to 3 imply supported, clamped, and 
elastically restrained edges are considered. 


INTRODUCTION 


The advent of high-speed flight and the concurrent necessity of 
maintaining aerodynamically smooth surfaces under high stress have led 
to the increased study of sandwich plate construction as a possible 
substitute for sheet— stringer construction in airplane design. A sand- 
wich plate consists essentially of a relatively thick, low-density, low- 
stiffness core bonded between two thin sheets of high— stiffness material. 
Materials that have been considered for the core include balsa wood, 
hard foam rubber, cellulose acetate, resin— impregnated cloth fashioned 
into a honey comb, corrugated metal sheet, and even closely spaced stiffeners 
of the conventional type. The face sheets may be of metal, plywood, 
wood-pulp plastic, or some other type of high— stiffness material. 

Because of the low— stiffness core, the sandwich plate will, in 
general, experience appreciable deflection due to shear. Furthermore, 
because the face sheets or core (or both) may have orthotropic stretching 
properties, the sandwich plate will in general be orthotropic in its 
flexural properties. As a result, ordinary plate theory, which is based 
on the assumptions that the plate is isotropic and that deflections due to 
shear may be neglected, cannot be used to determine the stresses, deflec- 
tions, or buckling loads of sandwich plates. 

A general small-deflection theory for flat orthotropic plates is 
therefore developed in which deflections due to shear are taken into 
account. The theory is applicable to any type of orthotropic or isotropic 
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sandwich that behaves essentially as a plate* provided certain physical 
constants are known. These physical constants (two flexural stiffnesses* 
two shear stiffnesses* a twisting stiffness* and two Poisson ratios 
defined in terms of curvatures) serve to describe the plate deformations 
associated with simple loading conditions and may be regarded as fundamental 
properties of the plate. For simpler types of sandwich construction the 
physical constants can be evaluated theoretically from the geometry and 
physical properties of the materials used. For more complicated types 
of construction* these constants can be evaluated by means of simple tests 
on samples of the assembled sandwich* as described in appendix A. A 
reciprocal relationship between the flexural stiffnesses and Poisson 
ratios is derived in appendix B. 

As is the case with ordinary plate theory* the orthotropic plate 
theory consists of two parts* each complete in itself. These parts are 
a set of six differential equations* three of which express the equilibrium 
of an infinitesimal plate element and three of which relate the curvatures 
and twist of the element to the forces and moments acting upon it* and 
an expression for the total potential energy of the system comprising the 
plate and the forces acting upon it. The six differential equations 
involve six variables. However* it is shown how these simultaneous 
equations can be reduced to a single equation of sixth order involving 
any one of the variables alone. In appendix C the consistency between the 
differential equations and the potential energy expression is shown by 
a variational method. 

The consideration of deflections due to shear makes necessary the 
specification of one more boundary condition than in ordinary plate 
theory. This fact was first appreciated by Reissner in reference 1. 

Because of some arbitrariness in the choice of the additional boundary 
condition* two types of simple support and two types of clamped edges 
are possible. Furthermore* three boundary conditions can be specified 
for a free edge* in contrast to ordinary plate theory. Boundary 
conditions more general than freedom* simple support* or clamping are 
considered in appendix C., 

A number of investigations related to the problem of orthotropic— 
or isotropic— sandwich— plate analysis have been made previously. Theories 
for the bending of orthotropic plates due to lateral loads and buckling 
due to edge loads* neglecting deflections due to shear* are given in 
references 2* 3* and- ^ and pages 38 O— 384 reference 5* The offset of 
shear on the bending of homogeneous isotropic plates and isotropic 
sandwich plates due to lateral load Is considered in reference 6 . The 
effect of shear on the bending due to uniform lateral load and buckling 
due to edge compression of simply— supported isotropic sandwich plates 
with homogeneous cores is considered in investigations by Hopkins and 
Pearson and by Leggett and Hopkins. A rough method of taking deflections 
due to shear into account In the buckling of simply-nsupported ortho- 
tropic sandwich plates is used in reference 7 - 
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Tlie present theory may be regarded as a natural extension to plates 
of the approximate theory used in pages 170 — 17 ^ of reference 8 to take 
into account deflections due to shear in a beam. The theory of this paper 
is more general than the aforementioned theories in that it applies to 
orthotropic or isotropic sandwich plates with homogeneous or nor>- 
homogeneous cores and with arbitrary boundary conditions, it presents 
both the differential ecuations and the energy expression for the plate, 
and it is applicable to problems that involve lateral as well as edge 
loads. The differential equations of the present theory are reduced to 
special forms in order that they may be compared with the equations 
obtained in references 5 and 6. 

The detailed development of the theory comprises most of the following 
sections and. the appendixes. The main parts of the theory are summarized 
briefly in a section entitled ,r I®CmTULATIOR OF PRINCIPAL RESULTS . " 


SYMBOL! 


y> z 


w 


<1 






^xy 


orthogonal coordinates; z measured normal to plane of 
plate and x and y parallel to principal axes of 

f lexural symne try , i nc lie s 

deflection of middle surface of plate, measured in 
z— di rec t i on, i nc he s 

intensity of lateral loading, pounds per square inch 

intensity of internal shear acting in z— direction in a 
cross section originally parallel to yz— plane, 
pounds per inch 

intensity of internal shear acting in z— direction in a 
cross section originally parallel to xz— plane, 
pounds per inch 

intensity of internal bending moment acting upon a cross 
section originally parallel to yz— plane, inch— pounds 
per inch 

intensity of internal bending moment acting upon a cross 
section originally parallel to xz— plane, inch— pounds 
per inch 

intensity of internal twisting moment acting in a cross 
section originally parallel to yz— plane or xz— plane, 
inch— pounds per inch 
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Nx intensity of middle-plane tensile force parallel to 

x/,— plane , pounds per inch 

Hy Intensity of middle -plane tensile force parallel to 

yz— plane, pounds per inch 


N- 


'xy 


Djj, Dy 


xy 


D 


intensity of middle-plane shearing force parallel to 
yz— plane and xz— plane, pounds per Inch 

flexural stiffnesses of plate with anticlastic bending 

unrestrained, inch— pounds ( 

y Curvature J 

twisting stiffness of plate, inch-pounds 

( Twisting moment /inch ] 

Twist J 

flexural stiffness of ordinary plate, inch-pounds 
shear stiffnesses of plate, pounds per inch 


hx> Py 


Poisson ratios for plate, defined in terns of curvatures 


P 


Poisson ratio for ordinary plate 


7 X , 7y 


shear- strain angles due to shears Q x and Qy, 
respectively, radians 


h 


thickness of plate, inches 


a, b 


length and width, respectively, of rectangular plate, 
inches 


V 

VI 

v 2 


total potential energy of system, inch-pounds 
strain energy of tending of plate, inch-pounds 
potential energy of external loads, inch-pounds 


u, v displacements in x— direction and y— direction, respectively 

of a point in middle surface of plate, inches 

[D], [M], differential operators 

[N], [P) 


SIGH CONVENTION 


The sign convention and notation used in the present paper Is, 
wherever convenient, the same as that used by Timoshenko in reference 5 . 
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The x—, y— , and z-axes of an orthogonal coordinate system are oriented 
so that the xy— plane coincides with the undistorted middle plane of the 
plate. Deflections w are measured normal to the xy— plane and are 
positive in the positive direction of the z— axis. The lateral load q 
is also positive in the direction of the z— axis. 

The internal shears and Qy, moments My* and “xy' 

and middle- plane forces N x , Ny, and ere shown in figure 1 acting 

in their positive directions upon an Infinitesimal element of length dx 
a nd width dy cut from the unloaded plate by planes parallel to the xz— 
and yz— planes. Only the forces and moments acting on two adjacent faces 
of the element are shown. The forces and moments on the opposite faces 
differ frcm those on the faces shown only by infinitesimal amounts . The 
directions in which they act, however, are opposite (e.g., mcment dy 

on the face shown is counterclockwise; moment dy on the opposite 

face would be shown acting clockwise). The twisting mcment and middle- 
plane shearing force acting on any cross section are known, frcm 
equilibrium considerations, to be equal to the twisting moment and middle- 
plane shearing force acting on a cross section at right angles. The 
symbols and therefore appear in both of the faces shown in 

figure 1. 

For convenience, in this report the z-direction is sometimes referred 
to as the vertical direction and planes parallel to the xy-plane are some- 
times referred to a3 horizontal planes. 


PHYSICAL CONSTANTS 


The physical properties of the plate are described by means of 
seven constants: the flexural stiffnesses D x and Dy, the twisting 

stiffness D^y, the transverse shear stiffnesses Dq^ and Pq^, and - 

the Poisson ratios p x and My. Definitions of these constants are 

obtained by considering the distortions of the differential element of 
figure 1 under simple loading conditions. 

Let all forces and moments acuing on the element be zero, except 
for the moments Mx acting on two opposite faces. Hie effect of is 

to produce a primary curvature in the middle surface of the element 

&x 2 

and also a secondary curvature which i3 a Poisson effect. Then D x i 

Sy^ 

defined a 3 the negative of the ratio of moment to primary curvature cr 


D-u 


dx 2 


(la) 
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when only is acting, and ix^ is defined as the negative of the ratio 

of Poisson curvature to primary curvature or 


he = " 


d^ 

d * 2 


(lb) 


when only is acting. No other distortions are assumed but 2 l3£ 

cDc 2 

and *-5 when he acts. The minu 3 signs are introduced in order to 

w 

make D x and ^ essentially positive quantities. 

Similarly, Dy and uy are defined as 


D y - 


& 

Sy 2 


(lc) 


hr 


& 

dx^ 

^w 

dy 2 


(Id) 


when only My is acting. 

If now all of the forces and moments are equal to zero except 
M^ pctlng on all four faces , the only distortion produced is a twist 

b 2 v 

— and Dw is defined as the ratio of twisting moment to twist or 

dx dy ** 


= -^ 2 - ( 1 .) 

d^w 
dx dy 


when only is acting. 
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The transverse shear stiffness Dq_^ la defined hy letting only 

the shears Q x act on opposite faces of the element (except for an 
infinitesimal moment of magnitude Q x dy dx required for equilibrium). 
The distortion is assumed for the moment to be essentially a sliding of 
one face of the element with respect to the opposite face, both faces 
remaining plane. As a result of this sliding, the two faces parallel 
to the xz— plane are distorted from their rectangular shape 
parallelograms by an amount y x , which is the Bhear angle 

the xz— plane. The shear stiffness Dq^ is defined as the 
shear to shear angle or 

r x 

when only Q x is acting. If the sides of the element are 
to the z— axis, the slope of the middle surface is 

dw _ _ Qx 

dx x Dq x 

when only Q^. is acting. 

In a similar manner, the shear stiffness Dq is defined as the 

ratio of the shear on the faces parallel to the xz— plane to the shear 
angle measured in the yz— plane when only Qy Is acting or 

<*> 

7 y 

wlen only Qy is acting. If all sides of the element are kept parallel 
to the z— axis, the slope produced is 

= 7 = 

dy y DQy 

when only Qy is acting. 

The constants Just discussed serve to define the orthotropic 
sandwich plate; they can be evaluated theoretically if the properties 
of the component parts of the sandwich are known and if the 
plate is of simple enough construction. In any event, the constants 
can be determined experimentally by means of bending tests and twisting 
tests on beans and panels of the same sandwich construction as the plate. 
A description of the tests required is given in appendix A. 

Although seven physical constants have been discussed, they need 
not all be Independently determined for if any three of the four 


into 

measured in 
ratio of 

(If) 

kept parallel 
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constants D^ Dy, \Xyr_, and Uy are known the fourth can be evaluated 
from the relationship 


This relationship, based on a generalization of Maxwell* s reciprocal law, 
is derived in appendix B. 


The shear stiffnesses Dq^ and Dq merit same additional discussion. 

The distortion due to shear was assumed to be a sliding of the cross 
sections over each other, the cross sections remaining plane and the shear 
strains remaining constant for the entire thickness of the plate and equal 
to the shear angle 7 X or 7 . Actually, if the plate is continuous enough 
for cross sections to exist at all, under shear the cross sections generally 
tend to warp out of their plane conditions (p. 170 of reference 8)j this 
warping makes the shear angle, as defined for equations (if) and (lg), 
meaningless. The shear strain varies with depth and an average shear strain 
will have to be used as the effective shear angle y x or 7 y for purposes 
of defining effective shear stiffness or Dq^. Id the experimental 

method is used (see appendix A), this difficulty is not encountered because, 
instead of a shear angle, curvatures are measured, and the stiffnesses 
obtained are automatically the effective stiffnesses. 


Despite the general tendency of cross sections under shear to warp, 
the assumption that they remain plane (though not normal to the middle 
surface) can be shown to be practically correct for those sandwiches in 
which the stiffness of the core is very small compared with the stiffness 
of the faces (e.g., metalite, honeycomb). For such sandwiches the shear 
stiffnesses Dq^ and 

may be assumed to take all the direct bending stress and the vertical shear 
may therefore be assumed uniformly distributed in the core. The shear 
angles y x and 7y will then be constant throughout the core. 




can be readily calculated, because the faces 


For those sandwiches in which cross sections under shear may not 
be assumed to remain plane, the tendency of these cross sections to warp 
introduces a further complication which can, however, be resolved by means 
of a Justifiable simplifying assumption. This complication is due to the 
fact that if the cross-sectional warping is partially or completely 
prevented the effect will be to increase the shear stiffness Dq x or 
Dq . The shear stiffnesses, thus, depend not only on the properties of 


the plate materials but also on the degree of restraint against cross- 
sectional warping. For the purpose of the present theory the 3hear 
stiffnesses Dq^ and are assumed to be constant throughout the plate 

and have the values they would have if cross sections were allowed to warp 
freely. The error caused by this assumption will be mainly local in 
character, being most pronounced in the region of a concentrated lateral 
load, whe re a sudden change in the shear tends to produce a sudden change 
in the degree of warping which is prevented by continuity of the plate. 

The error will probably be negligible in the case of distributed loads. 
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for which there are only gradual changes in the shear. A discussion of 
this error in connection with beams is contained in pages 173—17^ of 
reference 8 and in reference 9 . 


DIFFERENTIAL EQUATIONS FOR PLATE 
Distortion Equations 

^ 2 V A 2 V 

Equations can be derived relating the curvatures and — -5 

and the twist at any point in the plate to the internal shears 

dx dy 

and moments acting at that point. 

Equation for the curvature An expression can be obtained 

9 k 2 


for the total curvature 


2 bc 2 


in the x-dlrection by adding together the 


contributions made by each of the shears and moments acting separately: 
From equation (la) the curvature contributed by is found to be 


5k 

Dr 


d 2 w 


Equations (lc) and (id.) can be solved for the contribution to 
by 1 ^. which is 

Finally , "the equation following equation (if) indicates that the 
existence of 


_4-A> 


produces a curvature in the middle plane equal to 

1 dQn 
* ^x 


5^w 

dx 2 * 


The moment ; and the shear Qy make no contribution to 

Addition of tne three component curvatures gives 

dx 2 D x ^Dy D Qi dx 

Equation for the curvature ~ -5.— Similar considerations give 

oy d 


the curvature in the y-direction as 
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Sy 2 



Mx 

u x 


1 

’C 



(3b) 


Equation for the twist 


d^ 

bx by 


.— An expression for the twist 


b 2 * 

Sx~5y 

is obtained by first writing an expression for the twisting moment 
in terms of the distortions of the element dx dy. 


Let the middle surface of the element be distorted so that it 
acquires a twist J&L. Further assume that each line element normal 

&c dy 

to this middle surface before its distortion (a) first rotates so as 
to remain normal to it after its distortion, (b) then rotates through 
an angle y x in a plane parallel to the xz— plane, and (c) then rotates 

through an angle ?y in a plane parallel to the yz— plane . (Rotations 
(b) and (c) produce parallelogram-type distortions of cross sections 
and are therefore denoted as shear angles y x and 7y.) 


Distortion of the element as a result of rotations of type (a) is 
shown in figure 2(a). Distortion of the element due to rotations of 
type (b) is shown in figure 2(b) on the assumption that y x is zero 

at the center of the element and 1 b changing uniformly in the y— direction. 
Distortion of the element due to rotation of type (c) is shown in 
figure 2(c) on the assumption that y y is zero at the center of the 
element and changing uniformly in the x— direction. The magnitudes of 
the displacements shown on figure 2 are obtained "by considerations of 
geometry, the details of which are not given. 


The twisting moment M^y acting on all four cross sftstions of 
the differential element is proportional to the shear strain of the 
upper and lower surfaces, because this type of strain, JV "oughout the 
thickness of the element, produces the horizontal sb ^ couples 
that make up M—.. By superposition of the three d-.. s shown 


in figure 2, the shear strain in the upper (or lower; 
written as 


ce can be 


h d^ + h d 2 w \ 
2 dx by 2 bx by J 


h d?x __ h dy y 
2 "Sy ~ 2 bx 


or 


d^w _ 1 _ 1 

dx dy 2 dy 2 bx J 
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and therefore. 


m*, . h / -2%- - i |Zi - i £x 

^ \bx by 2 by 2 dx 


Substitution for y x and 7 y in terms of and Qy (equations (If) 

and (lg)) gives 


M-y ■ h* f • -S%- 

** \8* Sy 2 Hj,. W 2 Da ary 


where h' is a proportionality constant absorbing h. When Q x 
end Qy are both set equal to zero, the above equation must reduce 
to equation (le), because only Mxy is acting on the differential 
element. The constant h' is therefore identified as D^y and the 
equation for twisting moment becomes 


%y 



1 1 ^ 11 doA 

dxdy 2 Dq dy 2 Dn dx J 

JL *7 * 


Solution for 




yields the following equation analogous to the 
\o no 

equations already obtained for 212 and. 2_£; 

dx 2 dy 2 


J&. = ^ + + 

ax ay Dxy 2 Dq x by 2 Dq^ ax 


Equilibrium Equations 

The element dx dy must be in equilibrium under all the forces 
and moments acting upon it. This condition implies that certain 
relationships must exist among these forces and ncmente. These 
relationships can be derived by considering the changes that occur in 
the forces and moments frcm one face to the opposite and writing the 
equations of equilibrium for the element. The equations are the same 
as in ordinary plate theory. For equilibrium of forces in the 
x— and y— directions, these equations are obtained from equations (196) 
of reference 5 : 
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dx 


SNxy 

+ z-JSL - o 

by 


(Ua) 


— *■ + —22. = 0 (U 

by dx 

The equation for equilibrium of vertical forces is given at the top 
of page 305 of reference 5 as 



And the equation for equilibrium of moments about the y— and x-axes 
is obtained from equations (188) and (189) of reference 5 as 




by bx. 


^ bx by 


( 5 b) 

( 5 c) 


(Equations (5) are also derived in appendix C by minimization of the 
potential energy ») Note that the left-hand side of equation (5a) con, 
by virtue of equations (5b) and (5c), be simplified to 


52s + 5i 

ax »y 


If, as is customary in small-deflection theoiy, the middle-plane 
stresses V and N xy are assumed to be unchanged in the course of 

the plate’s deflection and equal to their initial values before 
application of lateral load, then equations (t) are automatically 
satisfied and equations (3) and (5) constitute the six fundamental 
differential equations that determine the forces, moments, and 
distortions throughout the orthotropic plate. They can be used in 
their present form or in the alternate form obtained in the following 
section. 
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Alternate Form of the Differential Equations 

The fundamental differential equations (3) and. (5) can he 
transformed so as to separate variables. Equations (3) are first 
solved for Mj, and M^y to obtain 


Mx = - 



b_ 

dx 


dv _ 



bw _ 'N 


(6a) 


My = - 




(6b) 


m™ = - 3^ I (6c ) 

^ 2** &\by \J dy\SE d Qx/J 

With the left-hand side of equation ( 5 a) simplified to and 

the above expressions for M^, and Mjy substituted into 

equations (5b) and (50), equations (5) beaone, after sane regrouping 
of terms. 





bx. by. 


w + 


(IeH 


(10^ = 



“ — * 

dx by 1 - d x dy 



+ 


1 5 m. d g + P*. ^5 — 1 

2 Sy2 (1 - h X hy )I ) Q X &C 2 


Ox 


1 . 
2 Dq^ 5x"5y 


^Xd y 

(! - P x Uy)I>Q y 5x dyj 


= 0 


+ 
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“ ^xy 


D„ 


dx 2 ay i - iixMy V x d* 2 9yVJ 


d3 6l 


i^:_a 2 


+ p y !i x 


|_ 2 Dq x Sx ^ (1 - MacMy;^ dx dyj 


Qx 


i3bc. H + 


a±__ 


D Qy dx 2 l 1 ,fi x u y^ D Oy 6y 2 


Qy = 0 


These three equations can he solved to obtain a differential equation 
for w alone in terms of q, an equation for Q x alone in terms of q, 
and an equation for Qy alone in terms of q. This separation is 
accomplished most easily, for the case in which N x , Ny, and N^y are 
constant throughout the plate, by treating the three differential equations 
as though they were algebraic equations and solving for w, Qx, and 
Qy by means of determinants. The terms in the determinants are the 
differential— operator coefficients of w, Q x , and Qy, appearing in 
the three equations. In expanding these determinants, the rule for 
multiplication of linear operators must be used. For example. 


dy 2 Sx dy dx dy3 

As a result of such a solution, the following differential equations are 
obtained for w, Q^, and Qyt 


[D]w 

= - [M]q 

(7a) 

[D]Qx 

= - [lllq 

(7b) 


= - [p]q 

(7c) 


where [D], [MJ, [N], and [P] are differential operators defined 
as 
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(continued on next page) 
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& 



H| OJ + 

tl 

p7 


With these substitutions made, equation (7&) becomes, after some transposition of tenns 
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d**w _ d^w d^+w 

r + 2 1 ■ —r + p 

dx 2 dy 2 dy^ 



^ 2 W d 2 w 

^ + N X 2 -^ + N V 2 -? +2N 3 

dx 2 Y dy 2 


\ 

dx by/ 


which is the same as equation (19T) of reference 5. 


■ The 

reference 6 

by use of Castigllano’s theorem of least work for the case in which the 
middle— surface forces N x , Ny, and N^y axe zero. The equilibrium 
differential equations of reference 6 are equivalent to equations (5) of 
the present paper. Equations ( 10 a), (lOd), (lOe), and (lOf) of reference 6 
can be solved simultaneously to obtain the following equations fox the 
curvatures and twist in terms of the vertical shears and moments (the 
notation is that of reference 6 ) : 


Isotropic sandwich plates, deflections due to shear considered 
di fferential equations for isotropic sandwich plates are obtained in 


, v% , 1 ^1 1 , av y \ 

ax 2 D(l-V 2 ) D(l-v 2 ) 0 8 dx C n \Sx Zy ) 



The symbols H, 
— , Qq£ f and 
The quantities 


V x , and Vy in the above equations correspond to 
Qy, respectively, in the notation of the present paper. 

D, C s , C n , and v are physical constants for the plate. 


The above equations are seen to be identical in form to equations (3) of 


the present paper (if Dq^. 
and y - directions ) except 


is set equal to D Qv 

y 

for the additional term 


for Isotropy In the x 

1 / av T 8v v \ . 


each curvature equation. This term arises frcaa the consideration of 
stresses and strains in the vertical direction, which were neglected in 


the present paper on the ground that they have a negligible effect on the 
over-all flexural behavior of the plate and are only important in the 


neighborhood of concentrated loads. Setting C n equal to infinity makes 
the equations derived from reference 6 completely identical in form to 
equations (3) of the present paper. It Bhould be mentioned that the 
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quantity w as used in reference 6 is not the deflection of the middle 
surface but "a weighted average across the thickness of the deflections 
of all points of the plate which lie on a normal to the middle surface." 


BOUNBART CONDITIONS 


The boundary conditions are first discussed for those types of 
edge support most commonly assumed in practice, namely complete freedom, 
simple support, and clamping. (More general kinds of support are 
considered in appendix C.) These supports are characterized by the 
condition that no work is done by the moments and vertical forces at 
the boundary. A boundary parallel to the y— axiB is considered; the 
conditions for a boundary parallel to the x-axie can be obtained by 
replacing x by y and vice versa, except in the subscripts of 
Mxy ^xy • 

Free edge .— The boundary conditions for a free unloaded edge 
parallel to the y-axis express the conditions of zero bending moments, 
zero twisting moment, and zero vertical force, or 


M* = 0 (9a) 

= 0 (9b) 

Qx = 0 (9c) 

If the free edge carries load, the middle— plane forces N x and 
will not in general be zero and the boundary condition of zero net 
vertical force becomes 


Qx + 


N — 

x Sx 


+ N. 


xy 


dw 

dy 


= o 


(9c r ) 


instead of equation (9c). 

Simply supported edge .— The principal boundary conditions for a 
simply supported edge parallel to the y— axis are w = 0 and Mx = 0. 

If to these two conditions is added the restriction that there is no 
y— displacement of points in the boundary, then the shear angle ?y is 

Qv 

zero and therefore =0. If, on the other hand, the support at the 

“Sy 

“boundary is applied only to the middle surface at the boundary and no 
horizontal forces are applied to prevent the y— displacement of other 
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points in the boundary, then M^, which is made up of such horizontal 
forces, must be zero. Two different types of simple support thus emerge. 


For simple support in which 
from moving parallel to the 

all points in the boundary are prevented 
edge, the conditions are 


o 

II 

> 

(10a) 


* 

II 

o 

(10b) 


o 

11 

"IS' 

(10c) 

For simple support in which all points 
in the middle surface, are free to move 
conditions are 

in the boundary, except those 
i parallel to the edge, the 


< 

II 

O 

(lla) 


- 0 

dib) 


Mxy = 0 

(lie) 


Of the two types of simple support, the first (equations (10)) is more 
likely to occur in practice. 

Clamped edge.— The principal conditions characterizing a clamped 
edge parallel to iiie y— axis are zero deflection of the middle surface 
and zero rotation of the cross sections making up the boundary 
(i.e., the boundary plane remains parallel to the z-axis) . The require- 
ment of zero deflection is satisfied by letting v = 0 at the boundaiy. 
The requirement that boundary cross sections remain parallel to the 

dw Q-v 


z— axis is satisfied by letting 


d q 


as the equation following 


equation (if) indicates, 
neglected by letting Dq 

.3w 


(Note that *.f deflections due to shear are 
= co, then the last boundary condition 


reduces to 




= 0, which is familiar in ordinary plate theory.) Just 


as in the case of simple support, the third boundary condition is 


or “zy = 0 depending on whether or not pointB in the 


Qv 

either — - 

% 

boundary (other than those points in the middle surface) are prevented 
from moving parallel to the edge. Thus, two types of clamping are 
possible. For a clamped edge in which the points in the boundary of 
the plate are prevented from moving parallel to the edge, the 
conditions are 
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v = 0 (12a) 

= 0 (12b) 

S D ^x 


Jfc - 0 (12c) 

D 0y 


For a cleuaped edge in which the points in the boundary (except those in 
the middle surface) are free to move parallel to the edge, the 
conditions are 


dx 


w = 0 

( 13 a) 

-Ss-0 

( 13 b) 



o 

II 

( 13 c) 


The latter type of damping is very unlikely to occur in practice, 
because any practical type of restraint that keeps the boundary from 
rotating has to be applied over an appreciable part of the thickness 
of the edge and therefore prevents most points in the boundary from 
moving freely parallel to the edge. 

The boundary conditions Just discussed, as well as boundary 
conditions corresponding to more general types of support, are 
derived in appe ndix C by a variational method. 


POTENTIAL- ENERGY EXPRESSION 
Strain Energy 


An expression can be obtained for the strain energy Vj produced 
by the moments M*, My, and M^y and the shears Q* and Qy by 
considering the work done by these moments and shears in distorting the 
differential element of figure 1. 

The work of the moments M x dy is equal to dy times the 

counterclockwise rotation of the right-hand face with respect to the 
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left-hand, face of the element. This rotation ia made up of two parts: 
the rotation caused by the moment Mj itself and the Poisson rotation 

caused by the moment My. The sum of these two parts is 

— ( — ^ (Note that although the term makes a contri— 

\ % •'Ey / dx 

bution to the curvature of the middle surface, this texm represents a 
rate of change of sliding rather than a rate of change of rotation and 
therefore makes no contribution to the rotation of one face with respect 
to the opposite.) The work of the moments M x 1 b therefore 


_ ^ dy [ - ^ + uy^ 




I dx 


or 


i: -JS& 

2x y iV 


) 


dx dy 


(a) 


Similarly, the work of the moments My. is 

i(^ - M ^) axaj 


(b) 


The work of those moments Mj-y acting in the faces parallel to 

to the xz— plane is equal to lx times the clockwise rotation of 

the nearer face (as seen in fig. 1) with respect to the farther face. 
This rotation is made up of the two parts shown in figures 2(a) and 2(b) 
and is equal to 


4 d7x j 

or, replacing y x by its equivalent in terns of Q x (equation (if)). 


c&r 
dx dy 


_ - \ 

d Qx / 




The work of the moments parallel to the xz— plane is therefore 
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1 


D Qx 



dx dy 


Similarly, the work of those moments V parallel to the yz— plane 1 b 



5^w 
bx. by 


-i- ay 

\ * J 


Tlie total work of the moments M^y is, by adding the last two 
expressions , 


Mxy' 


& _ i -L ^ -1 -L ay 

chTciy 2 Dq x by 2 D Qy c»x J 


The factor in paxentheses is simply , from the equation preceding 

u xy 

equation (3c), and the work of the moments M^y therefore becomes 



dx dy 


(c) 


The work of the shears Q x is ^Qx dy times the downward 

distance through which the right-hand fane slides with respect to the 
left-hand face. This distance is ? x dx and work is therefore 

■1q,x7x dx dy . Replacement of y by its equivalent in terns of Q x 
2 * 
gives 

, n 2 

— -2— dx dy (d) 

2 \ 


for the work of the shears Q^. Similarly, the work of the shears Qy 
is 


1 % 

2 £>Q 


dx dy 


V 


!e) 


Integration of the energy expressions (a) to (e) over the entire 
plate gives, as the total strain energy due to bending and shear. 
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The potential energy acquired by the external forces in the course of the lateral deflection 
of the plate is independent of the internal construction details of the plate and depends only on 
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"by addition of equations (lM) and (15)> 
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d 2 * 


d 2 ^ 

dy 2 


a 2 * 

cir dy 


= _*k 


Dr 


+ Uy 


% 


+ 


_L?2i 

ax 


= _!k 
% 


+ U* 


*k 


DQy 


= «c + i±^ + ii.^ 

V 2 Dq x dy 2 DQy dr 


relating distortions to distorting moments and forces, and 



d 2 v 

1 + ^2 + N y 


d 2 w d 2 v \ 

^ + 2 K ^aTd 7 y 


Ox 


*Jki , h Jk 

dy dx 


a _a ^ 2 + a^ 

. ax a y 


for equilibrium. 

3. The first three equations can be solved for 
May to obtain 


Mr* 


V 


and 


Mr -- 


Dr 

1 ~ ^^y 


My =- 




Mxy ■ |Dxy 


dr 


M _ Bi_\ + d_/ dst _ 
dy Dq^ J ^y\dr Dq x 



Substitution of these expressions into the last three equations and 
solution of the resulting equations by means of operational determinants 
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gives the following differential equations with variables separated, for 
the case in which N XJ Ny, and N^y are constant throughout the plate: 

[D]w = - [M]q 

[D]Qr = - [N]q 

[D]Qy = - [p]q 

where [D], [M], [N], and [P] are differential operators defined 

by equations (8) , 

4. Three types of support commonly assumed at the boundaries of 
a plate are no support (free edge), simple support, and clamping. 

These types of support can be described in terms of deflection, shears, 
and moments for an edge parallel to the y— axis as follows: 

For a free edge. 


M* = 0 


Q, 


M xy 

N 

x 


= 0 


+ N 




8w 


= o 


For a simply supported edge at which the support is applied over the 
entire thicknesB, 


w = 0 
= 0 



0 


For a simply supported edge at which the support is applied only to 
the middle surface. 


w = 0 
M x = 0 


“xy = 0 
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For a clamped edge at which the support is applied over the entire 
thickness. 


w = 0 


^w Qx _ Q 

Si Dq x ' ~ 

% 


For a clamped edge at which the support is applied only to the 
middle surface (a type of support very unlikely to be met in practice). 


w = 0 

aw _ Qx_ = 0 
c* »Q X 

”xy - 0 


The conditions for an edge parallel to the x— axis can he written by 
replacing x by y and vice versa, except in the subscripts of 

^xy ^xy • 


Boundary conditions can also be written for mol's general types 
of support. (See appendix C,) 

5 . The potential energy of a plate in which the middle-surface 
forces are assumed to remain unchanged in the course of the plate f s 
deflection and for which the moments and vertical forces at the 
boundaries do no work is 
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APPENDIX B 

DERIVATION OF RELATIONSHIP d x D y = Myl> x 


Betti’s reciprocal theorem (reference 10) can he expressed as 
follows: Let two groups of forces he applied to a structure, each 

group of forces producing distortions that are directly proportional 
to the magnitude of the forces; then the work of the first group of 
forces anting through the displacements produced hy the second group 
is equal to the work of the second group acting through the displace- 
ments produced hy the first. 

The structure to which this principle is applied is the element 
cbc dy of figure 1. Let the first group of forces consist of the 
moments M^- dy. The distortions produced are the curvatures 

and where, from equations (la) and (ih), 

dx 2 dy 2 


and 


5 2 v _ Mx 
Sx 2 D x 


5y 2 



The second group of forces are the moments 
produced by the group are the curvatures 
equations (lc) and. (id). 


My djc, 

£2 and 

dy 2 


and the distortions 

- — where, from 

dx 2 


and 


ctw _ 

*y 2 D y 





£ 


The work done hy the first group of forces M x dy in association with 
the curvature n v £ produced hy the second group is 

* T\ 



3^ 
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or 



dy 


Similarly, the work done by the second 


association with the curvature 


group of forces My dr In 
produced by the first group Is 


-My dx^gS dy) 
or 


- ! VVn" ** ** 

v x 


Equating the expressions for the two works and eliminating the common 
factor — M^My dx dy gives 


iix = iis 


from which is obtained equation (2). 
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Minimization of total potential energy .- The conditions that must be satisfied if the total 
potential energy V of the system is to be a minimum are now considered. By the calculus of 
variations (reference 11 ) , minimization of V requires the vanishing of the first variation 57 
The first variation can be evaluated from equation (Cl) as 
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The first of these two integrals can be integrated again by parts so as to leave only terms 
containing 5w rather than derivatives of 5w. With the aforementioned integrations by parts 
performed and after considerable rearrangement of terms, the expression for 8V becomes 
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By virtue of equations (6) and ( 2 ), the above expression for 5 V can be revrritten as 
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differential equations result from equating the surface integrals to zero: 
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^ - 2 ¥0SL + ^ = -( q + N x ^ 
clx 6y 6y^ \ 6x‘ 


Sx 6y 


dy dx 

*“a + ^ 

Sx Sy 


By virtue of the last two of equations (C2), f ) and. 

V dx dy / 

( — ffixy J in the line integrals can "be replaced by Q™. and Qy, 

\oy ox J J 

respectively. Equating the line integrals to zero then gives the 
following boundary conditions required to insure that 5V = 0: 


Qx + Nx S + *** SF = ^ X ° r ^ = ° 


My = M x or &(§& - ) = 0 

\Sx 0 q x/ / 


1% =M ^ ^ \^~ 0 tj = ° 


°y + % ^ + N xy ^ = Q y 


or 6w - 0 


° r 


or 


at y = 0,b 
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Equations (C2) are the differential equations that must be 
satisfied if the potential energy is to be a minimum. They will be 
recognized as the equations of equilibrium, equations ( 5 ). 

Equations (C 3 ) and (C4) are the boundary conditions that must 
be satisfied if the potential energy is to be a minimum. The lefi>- 
hand groups of equations (C3) and (C4) imply that the limiting values 
of the internal forces and moments, as the edge of the plate is 
approached, must be in equilibrium with certain prescribed forces and 
moment s externally applied at the edge (the prescribed forces and 
moments being designated by means of the horizontal bars). The right- 
hand groups of equations (C 3 ) and (C4) imply that the displacements at 
the edge must have certain prescribed values. 


The boundary conditions given by equations ( 9 ) to ( 13 ) for free, 
simply supported, and clamped edges parallel to the y-axls are special 
cases of equations (C 3 ). For example, the boundary conditions for a 
simply supported edge (equations (lO))can be obtained from equation (C 3 ) 


by prescribing the values of w, Mr, 
at the boundaries x = 0 ,a. 



to be zero 


If a plate is elastically supported at the boundaries, the elastic 
support may sometimes be conveniently thought of as made up of three 
rows of closely spaced discrete springs at each edge: a row of 

deflectional springs, a row of rotational springs, and a row of torsional 
springs, having the known stiffnesses per inch kq, k 2 , and k^, 

which may vary along the edge. For this type of support the vertical 
shear reaction at any point along the edge is proportional to the 
vertical deflection at that point and the twisting and bending moment 
reactions are porportional to the corresponding rotations of an 
originally vertical line element in the edge. The boundary conditions 
for this type of support can be obtained frcm equations (C 3 ) and (C4) 
by setting 


Ox = kpw at x = 0 


§ x = — kjw at x = a * 



at 


x = 0 
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'■-‘‘(fc-Sfc) 

= kjw at y = 0 
Qy = - k^y at y = 

a. = - - 2z- 

^ V* % 

* " *(* " \) 

^(MO 


at x = a 


at x = 0 


at x = a 


: b 

^ at y = 0 

at y = b 
at y = 0 


^ at y = b 


The signs in the above boundaiy conditions follow as a result of the 
directions assumed for positive shears and moments. 
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APPENDIX D 

DERIVATION OF EQUATION (15) FOR THE POTENTIAL ENERGY 
OF THE EXTERNAL FORCES 


A rectangular plate the edges of which are x = 0,a and 
y = 0,b is considered (fig. 6). The boundary conditions assumed 
are the usual conditions corresponding to zero work by the reactions; 
that is, each edge Is either free, simply supported, or clamped. 

The horizontal loads N x , Ny, and Nxy are assumed first to 
be applied at the boundaries with no lateral load. As a result the 
middle plane (and all' horizontal planes) of the plate stretches; thus, 
the constant stretching energy discussed previously in connection 
with the strain energy of the plate is produced, and slight shifts 
in the points of application of the edge forces N x , Ny, and N^ 

are caused. These new positions of the points of application are 
used as the arbitrary fixed reference points in any future measure- 
ments of the potential energy of the horizontal edge forces. 

If the lateral load q is now applied, the middle surface acquires 
the displacements w(x,y) In the z— direction, u(x,y) in the 
x— direction, and v(x,y) in the y— direction. As a result of these 
displacements, the lateral load acquires the potential energy 



qw dx dy 


(a) 


and the edge forces acquire the potential energy 



E x u 


0 


Nxyv 




NxyU 


b 

0 


) 


dx 


0 >) 


The moments and vertical forces at the boundaries do no work and 
therefore acquire no potential energy during deflection. 

By use of the formula for integration by parts > expression (b) 
for the potential energy of the ed.ge forces can be rewritten in terms 
of the interior forces and displacements as 
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In the development of the differential equations and in this section 
the middle-surface stresses N x , N y , and N^y are assumed to remain 

unchanged in the course of the plate’s deflection. Equations (4) for 
equilibrium of horizontal force, consequently, remain satisfied at 
all times, and, therefore, the last two integrals of expression (c) 
vanish. Furthermore, the assumption that the middle— surface stresses 
remain unchanged implies that no stretching of the middle surface 
during deflection occurs. In order to prevent Buch stretching the 
horizontal displacements u and v can be shown (p. 31-3» reference 5 ) 
to be related to the vertical displacements w as follows. 


du 1 /dv \ 2 

35E =“ 2^ ) 


bv _ l/aw\ 
by 2 \ oy ) 

Ju b* b w 

^ + <5x = ~ 3x ciy 


The first and only remaining integral of expression (c) therefore 
becomes 


1 

2 



+ 2N. 


d w 


5x 



dx dy 


(d) 
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Addition of expressions(a) and (d) gives, as the total potential 
energy of the external forces. 


V 2=5 Jf <M) 


Although the derivation was carried out for the special case of 
a rectangular plate, equation (Dl) also applies to a plate of any 
shape in which the middle— surface stresses remain unchanged during 
deflection. Equation (Dl) is identical with equation (15). 
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(a) Distortion due (b) 

t0 sTaP a11 


line elements 
remaining normal 
to middle surface. 


Distortion due 
to - * * . 

dy 


(c) Distortion due 

d7 v 

to TT- 


Figure 2.- Distortions of element dx dy in twisting. 
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